The Full Mottness 



Tudor D. Stanescu 

Department of Physics and Astronomy, Rutgers University, 
136 Frelmghuysen Rd., Piscataway, NJ., 08854-8019 

Philip Phillips 

Loomis Laboratory of Physics, University of Illinois at Urbana-Champaign, 1100 W.Green St., Urbana, IL., 61801-3080 

Though most fermionic Mott insulators order at low temperatures, ordering is ancillary to their 
insulating behaviour. Our emphasis here is on disentangling ordering from the intrinsic strongly 
correlated physics of a doped half-filled band. To this end, we focus on the 2D Ifubbard model. 
Because the charge gap arises from on-site correlations, we have been refining the non-perturbative 
approach of Mancini and Matusumoto— which incorporates local physics. Crucial to this method is 
a self-consistent two-site dynamical cluster expansion which builds in the nearest-neighbour energy 
scale, J. At half-filling, we find that the spectral function possesses a gap of order U and is devoid of 
any coherent quasi-particle peaks although ordering or charge fractionalization are absent. At low 
temperatures, local antiferromagnetic correlations emerge. In the doped case, we find that the Fermi 
surface exceeds the Luttinger volume. The breakdown of Luttinger's theorem in the underdoped 
regime is traced both to the dynamically generated Mott gap as well to a non-vanishing of the 
imaginary part of the self energy at the Fermi level. Spectral weight transfer across the Mott gap 
also emerges as a ubiquitous feature of a doped Mott insulator and suggests that high and low 
energy scales are inseparable. Additionally in the underdoped regime, we find that a pseudogap 
exists in the single particle density of states as well as in the heat capacity. The pseudogap (which 
is set by the energy scale t^ /U) is argued to be a ubiquitous feature of a lightly doped Mott state 
and simply represents the fact that hole transport involves double occupancy. In analogy with the 
Mott gap and antiferromagnetism, we propose that ordering may also accompany the formation 
of a pseudogap. We suggest a current pattern within a 1-band model that preserves translational 
but breaks time-reversal symmetry along the canonical x and y axes but not along x — izy that is 
consistent with the experimental observations. Finally, we show that the Hall coefficient in a doped 
Mott insulator must change sign at a doping level x < 1/3. The sign change is tied to a termination 
of strong correlation physics in the doped Mott state. 



I. INTRODUCTION 



Electronic systems with an odd number of particles per unit cell are typically metallic at zero temperature. However, 
Mott^ proposed that such systems in the presence of strong interactions can insulate at zero temperature without any 
accompanying symmetry breaking process, such as antiferromagnetism or charge ordering, which necessarily double 
the unit cell. Strictly then, a Mott insulator (should it exist) is a paramagnetic state with an odd number of particles 
per unit cell. Insulating behavior arises from the charge gap generated by the projective mismatch between the sub- 
lattices which have zero or some finite fraction of doubly occupied sites. When the overlap between such sub-lattices is 
sufficiently small, no transport obtains. In bosonic systems with a single boson per site, a true zero temperature Mott 
insulating state is realized in the standard quantum rotor model or Bose-Hubbard model when the on-site charging 
energy exceeds a critical value such that phase coherence is destroyed^*^. In fact, the recent observation^ that a 
Bose condensate in an optical lattice can be tuned between a superfluid and a Mott insulator simply by changing the 
intensity of the laser light places the bosonic Mott insulator on firm experimental footing. 

However, for fermionic systems, the inherent problem with the Mott insulating state of a half-filled band in d > 1 is 
its proclivity to order at zero temperature. Consequently, it is tempting to equate the Mott insulator with the ordered 
state or to assert categorically that the Mott insulator as an entity distinct from a symmetry-broken state never existed 
in the first place. In the context of high Tc in the cuprates, which are all anti-ferromagnetic Mott insulators, both 
views^ have been strongly expressed. In fact, a large body of work on the cuprates has focused primarily on models 
that capture low-energy spin physic o'^i^i^ at or near an antiferromagnet or a charge-ordered stat9i2iii4iSiiiMiiS4iSiiLi&, 
or a classificationi^j^S of the various charge-ordered states that ensue in a half-filled band. Alternatively, low-energy 
spin liquid niodel32iiS242i2ii2SiSS42L2& (that is, models with spin translation and spin rotation symmetry) have been 
proposed as candidates to describe the Mott insulator. In such approaches, the high energy scale associated with the 
charge gap is argued to be irrelevant, hence the focus exclusively on the spin sector to characterise the Mott insulator. 
Of course such an approach presupposes that the high and low energy degrees of freedom can be disentangled. 
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Should the insulating state in a half-filled band prove to be nothing other than a mean-field broken-symmetry state, 
then fermionic Mott insulators do not exist. Hence, a relevant question for the Mott insulating state is as follows: If 
we subtract the fact that ordering obtains at zero temperature, is there anything left over that is not explained by 
ordering? Equivalently, does ordering provide an exhaustive explanation of the state proposed by Mott? We refer 
to whatever might be left over once we subtract the fact that ordering has occurred as Mottness. That something 
might be left over is immediately evident from the fact that charge and spin ordered states all result in a doubling 
of the unit cell and hence are adiabatically connected to an insulator with an even number of electrons per unit cell. 
A typical example of such a system is a band insulator. Contrastly, Mott insulators which have an odd number of 
electrons per unit cell are not. Additionally, spin and charge ordered states can be described at the level of Hartree- 
Fock simply by constructing the correct broken symmetry ansatz. However the Mott state has no weak coupling 
or Hartree-Fock counterpart. Hence, although the Mott state might be unstable to ordering at low temperature, 
some features should be left over which represent the fingerprint of the nonadiabaticity with a band insulator and 
the fact that this state arises fundamentally from strong electron correlation. For example, experimentally22i2£, it is 
clear that above any temperature associated with ordering in both the electron and hole-doped cuprates, a charge 
gap of order 2eV is present in the optical conductivity and oxygen K-edge photoemissionp^. Hence, the vanishing 
of the low-energy (< leV) spectral weight at high temperature is not linked to magnetism or ordering of any kind. 
Further, the electronic bands below and above the charge gap are not rigid as would be the case in a band insulator. 
To illustrate, as a function of doping, in both hole and electron-doped cuprates^^'^*^ , the low-energy spectral weight 
increases at the expense of the high-energy (> 2eV) spectral weight such that the total integrated optical conductivity 
remains constant up to AeV. The same massive reshuffling of spectral weight from 2eV above the Fermi energy is also 
observed in 1-particle probes such as oxygen K-edge photoemissionSi and angle-resolved photo-electron spectroscopy 
(ARPES)22iS. Such spectral weight transfer indicates that the total number of low-energy degrees of freedom in 
the normal state of the cuprates cannot be decoupled from the high energy scales. What is surprising about the 
cuprates is that even when superconductivity obtains, the low and high energy degrees of freedom are still coupled. 
For example, Riibhaussen, et. ali^ have shown that changes in the optical conductivity occur at energies 3eV (roughly 
lOOA, A the maximum superconducting gap) away from the Fermi energy at Tc, and Bolegraf, et. alr-^' have seen 
an acceleration in the depletion of the high energy spectral weight accompanied with a compensating increase in 
the low-energy spectral weight at and below the superconducting transition. Similarly, Bontcmps, et. al,"^^ have 
directly observed that in underdoped (but not overdoped) BSCO, the Glover- Ferrel-Tinkham sum rule is violated 
and the optical conductivity must be integrated to 20, OOOcm"^ to recover the spectral weight lost upon condensation 
into the superconducting state. In a standard BCS superconductor, condensation leads to loss of spectral weight 
at energy scales no more than ten times the pairing energy. The fact that pair condensation perturbs the optical 
conductivity on energy scales as large as lOOA suggests that there is a direct link between the high energy Mott scale 
and superconductivity. Further, the persistence of spectral weight transfer through the superconducting transition 
indicates that the properties of the Mott state remain intact even in the presence of ordering. While spin-density wave 
antiferromagnets also possess two bands with a gap, such a state is insufficient to explain the origin of the spectral 
weight transfer in the cuprates. The reason is simple: spectral weight transfer persists well above T/v and at a doping 
level ^ 2%, where antiferromagnetism is absent. 

Concretely, what the optical conductivity and oxygen k-edge photoemission experiments on the cuprates lay plain 
is that regardless of whether ordering obtains at sufficiently low temperatures, the charge degrees of freedom remain 
characterized by a distinctly different state than that of any ordered state or a band insulator. To understand what 
physics is entailed by this state, consider the Mott mechanism for generating an insulating state in the one-band 
Hubbard model. At half filling, the chemical potential lies in the middle of the gap separating the lower and upper 
Hubbard bands which are dynamically split by the on-site energy for double occupancy. However, half the spectral 
weight resides in each of the bands. Consequently, to satisfy the sum rule that each state in the first Brillouin Zone 
(FBZ) carries unit spectral weight, the spectral function must be integrated across the charge gap not simply up to 
the chemical potential. Hence, the half-filled state is characterized by the Fermi energy lying in a gap but partially 
occupied states exist. It is this seemingly contradictory state of affairs that is at the heart of Mottness. Spectral 
weight transfer cannot obtain without it. For example, if each state below the chemical potential had unit spectral 
weight, no state would be available for spectral weight transfer from high energies. As a consequence, adding or 
removing an electron cannot be done without affecting both high and low energy scales. Consequently, at any doping 
level, the electronic states describing the charge carriers can be written as linear combinations of excitations living 
in both the LHB and UHB as will be detailed below. As a result, in the Mott state, the traditional notion that 
the chemical potential demarcates the boundary between zero and unit occupancy fails fundamentally. Of course, in 
Fermi liquids, the spectral function for each k— state can also have an incoherent background which can extend to 
high energies. However, as long as a coherence peak exists, a sharp criterion exists for unit occupancy of each state, 
namely whether or not the coherence peak crosses the Fermi level. In a Mott insulator, no such coherence peak exists 
and consequently, incoherence dominates the Mott state. 
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Doped Mott Insulator 





FIG. 1: Spectral weight transfer in a doped Mott insulator. The Photoelectron spectrum (PES) denotes the electron removal 
states while the electron-addition states are located in the inverse photo-electron spectrum (IPES). The on-site charging energy 
is U . Removal of a single electron results in the creation of two single particle states at the top of the lower Hubbard band. 
By state conservation, one state comes from the lower and the other from the upper Hubbard band and hence spectral weight 
transfer across the Mott gap. 



Alternatively, one can view the spectral weight transfer in real space by simply counting the number of available 
states for the photo-emission and inverse photo-emission spectra as demonstrated by Meinders, Eskes, and SawatzkjiSi. 
We recount the argument here as it is simple and instructive. Consider the half-filled one-dimensional chain of one- 
electron atoms shown in Fig. (Q. Both the electron removal (photo-electron) and electron-addition (inverse photo- 
emission) spectral weight are equal to N because at half-filling there are N ways of adding or subtracting an electron 
from a site that is singly occupied. When a single hole is added, both the electron-removal and the electron-addition 
spectral decrease to — 1 as there are now — 1 ways to add or subtract an electron from sites that are already 
occupied. Hence, there are two less states. These states correspond to the spin up and spin down states of the empty 
site and hence belong to the LHB. Consequently, the low-energy spectral weight (LESW), A(a;), has increased by two 
states. One of the states must come from the UHB as the high energy part now has a spectral weight of A^ — 1. Hence, 
for a single hole, there is a net transfer of one state from high to low energy. This argument is simply the real-space 
restatement of the more general principle that in k— space the LHB and UHB arc not static but dynamic and hence 
necessarily give rise to spectral weight transfer. In general, simple state counting yields 2x for the growth of A(a;) 
and 1 — x for the depletion of the high energy sector. In actuality, the dynamical contribution to the LESW results 
in A(x) > 2x. The dynamical LESW corresponds to virtual excitations to the UHB. Hence, in a strongly correlated 
system, the phase space available to add a single electron exceeds the nominal number of states initially present in 
the low energy scale, leading thereby to an inseparability of the low and high energy scales. 

The ubiquity of spectral weight transfer in the cuprates places extreme restrictions on which low-energy models 
are valid. Consider the standard t — J model^^. As the t — J model projects onto the LHB, we can estimate the 
corresponding LESW exactly by counting the electron-removal states. Consequently, A{x) = 2x is exact in the t — J 
model. However, in the actual Hubbard model, A(x) > 2x. Consequently, the standard t — J model does not have the 
correct number of low-energy degrees of freedom to describe the low-energy physics. This problem can be corrected2&, 
however, by using the full strong-coupling Hamiltonian that results from the t'^/U expansion and replacing all the 
electron operators by their projected counterparts. However, the price one pays is that the new projected electron 
operators do not obey the standard fcrmion commutation relations. That a correct low-energy theory must give up 
on either standard fermion commutation relations or particle conservation has already been pointed out by Meinders, 
Eskes, and Sawatzky^i. Hence, all standard low-energy fermionic theories do not have the correct physics to describe 
the cuprates primarily because the low-energy degrees of freedom are not fermionic. 

Given the ubiquity of spectral weight transfer over a large part of the phase diagram of the cuprates, it is imperative 
that any theory of high Tc incorporate the high energy scale associated with the charge gap. Further, because ordering 
does not seem to be a requirement for spectral weight transfer, our focus is on an accurate description of the high- 
temperature charge vacuum of the doped Mott insulating state. It is crucial in such an approach that the hierarchy 
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of energy scales, U, t, t^/U, for example, emerge. Consequently, we have been refining the non-perturbative approach 
of Matsumoto and Mancini^ to describe the inter-relation between the energy scales in a doped Mott insulator. This 
approach is based on the experimental observation that the low energy scales in doped Mott insulators are derived 
from high energies. As a result, the beginning point, namely the Hubbard operators, is one that is well-known to 
yield the Mott charge gap scale. Successively smaller energy scales are derived by treating local correlations on a 
small cluster. The self-energy of the lattice is then determined self-consistently from the local impurity problem. The 
Matsumoto-Mancini'*^ approach is then in the spirit of the cellular dynamical mean field treatment We report here 
the full details of this approach and catalogue the general results that follow from Mottness. Aside from a LESW that 
exceeds the nominal value obtained from state counting, we find that 1) Mottness gives rise to broad spectral features 
in the underdoped regime, 2) a violation of Luttinger's theorem in the underdoped regime, 3) hole-like Fermi surface 
near half-filling, 5) a jump in the chemical potential upon doping and 6) a pseudogap in the underdoped regime 
without invoking any symmetry breaking. In general, we find that the pseudogap is due strictly to near-neighbour 
correlations and can be thought of as the nearest-neighbour analogue of the on-site generated Mott gap. 

II. METHOD 

Many years ago, Hubbard--'' wrote the electron annihilation operator, Ci„ as a linear combination 

Cicr = Via +^tcr (1) 

of two composite excitations that reflect the energetically challenged landscape an electron must traverse in the 
presence of a large on-site Coulomb repulsion, U. Physically, the operators rji^ = Ci^Ui-a- and ^io- = ^^(l ~ 
represent an electron moving on doubly and singly occupied sites, respectively. Because such sites are split by C/, the 
Hubbard operators lead naturally to a gap at half-filling in a paramagnetic state, a result which thus far, only two other 
methods, dynamical mean-field theories (DMFT)-ii and quantum Monte Carlo (QMC)^2il^, have been able to obtain. 
Regardless of this success, the Hubbard operators have been criticized extensively because untested approximations 
generally accompany their implementation, and they lead to a Fermi surface which does not preserve the Luttinger 
volume. However, methods that are deemed as reliable^^i^^i^S^i, also find that in the lightly doped regime, Luttinger's 
theorem is violated as is seen experimentally^. Hence, a violation of Luttinger's theorem is not an a priori reason to 
dismiss the Hubbard operators. In fact, any violation of Luttinger's theorem must occur in the lightly doped or non 
perturbative regime. The untested or uncontrolled approximations usually arise from truncations in the equations of 
motion. However, such problems can be circumvented by the following procedure. First, project all new operators 
that arise from the Heiscnberg equations of motion of the Hubbard operators onto the Hubbard basis. Second, write 
the self energy exactly in terms of the remaining operators which are now orthogonal to the Hubbard basis. Third, 
use local DMFT methods to calculate the resultant electron self energy. The approximation introduced in the third 
step is that the self energy for a finite cluster is used to determine the self-energy for the interacting lattice. However, 
such methods have been shown to be strongly convergent and in fact constitute the accepted methodology for treating 
strongly correlated systems. In principle, as the cluster size is extended to infinity, such a cluster procedure becomes 
exact. Hence, an implementation of the Hubbard operators coupled with DMFT-type technology places the limitations 
not on truncation in the equations of motion but on the accuracy of the impurity solver and the size of the finite 
cluster. It is such a procedure that we outline here. As many of the details have been left out in previous presentation 
in the original paper by Mancini and Matsumotc^^ on the two-site cluster and subsequent implementations^, we will 
provide a complete derivation of the method so that any one reading this paper can implement it immediately. The 
key features of this method are its ability to describe physics on the scale of the Mott gap, [/, as well as on the the 
scale J « t^/U. 

A. Dynamical Green Function 

Our starting point is the on-site Hubbard model 

with nearest-neighbor hopping, tij — taij. We also introduce the composite operator basis 



(2) 



(3) 
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and its associated retarded Green function S{i,j,t,t') = {{ilJicr]'4'jc)} = 0{t — t'){{ipia{t),ipj^{t')}) . Writing the 

equations of motion for the Hubbard basis and projecting with the Roth^*' projector 7^(0) = S/n({^' ''/^/})-^;ri^V''nj we 
obtain for the 'current' operator the expression 

j,(t) = i^^, = K^j, + V{5],) + SJ, = E^, + 5J,. (4) 
The formal sohition for the Green function in Fourier space, 

contains the self-energy E(k,a;) = (5m(k, cj)I^^ with 

6m{\i,uj) = FT{9{t^t'){5J{t),5j\t')})i (6) 

where the subscript / indicates the irreducible part and FT denotes the Fourier transform. We consider the param- 
agnetic case, for which the overlap matrix I is diagonal with In = Ii — 1 — n/2 and I22 = I2 = n/2. 

The primary operational hurdle is the evaluation of the dynamical correction 6m.. The first step is to write the 
explicit expressions for the operators SJ which are 'orthogonal' to the basis ip. Using the notation t = 2d ■ t, d being 
the dimensionality of the system, and p ^ p — /A, we obtain 

(z) = ~i [tt,, + '^ct - ~ - Pi^Ii' - C^2"')] (7) 

and SJ2{i) — —5Ji{i), where the self-consistent parameters e and p = p + as well as the higher order composite 
operator tt^ are given by 



P = (ni^<^) (4^Cix(cJ^c,|)") - (citCix(c|^cJ^)"), 

and the superscript a denotes the averaging over nearest-neighbor sites. Consequently, we can write the dynamical 
correction matrix (5m in the form 

(5m(k, Lo) = Dm(k, uj) ^ }^ ~^ ^ (9) 

and the problem reduces to the determination of the higher-order Green function Dm(k) — FT{9{t — 
t'){SJn{t)5Jli{t)}) I . Because £'TO(k, w) cannot be evaluated exactly, we seek a systematic way of calculating the 
dynamical corrections. The simplest approach would be to consider the single-site approximation. Such an approxi- 
mation is in the spirit to the methods, in which the self-energy is momentum independent. An improvement 
would be to consider the dynamics associated with two sites as proposed by Mancini and Matsumoto^. Evaluation 
of the self-energy over successively larger clusters would lead to an exact determination of the dynamical corrections. 
The essence of this approach is based on the fact that the physics of strongly correlated electrons emerges mainly 
from local correlations: on-site interactions generate the Mott gap, while nearest-neighbour interactions generate the 
scale, /U . e must treat all excitations on equal footingSuccessively larger clusters build in lower and lower energy 
scales. However, due to the fact that transfer of spectral weight from high to low energies is characteristic of Mott 
insulators, a separation of energy scales in strongly correlated problems is not possible. Therefore, it is crucial that all 
the excitations be treated on equal footing and, consequently, local correlations must be included. Following Mancini 
and Matsumoto^, we write the dynamical corrections as a series 

Dm{x, x') = Sx.x'DmQ{x, x') 4- ^ Sx+a,x'D'mi{x, x') -\ (10) 

a 

in increasing cluster size. Here, x and x' are two representative sites and a indexes all nearest-neighbor sites. In the 
two-site approximation, the series is truncated at the level of on-site, Dttiq, and nearest-neighbor, Dmi contributions. 
In Fourier space, the dynamical corrections can be written as, 

Dm{k,uj) K Dmo{uj) +a{k)Dmi{uj). (11) 
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Here Dmo and Dmi involve Green functions of operators defined on nearest neighbor sites (see Eq. [Tj). Consequently, 
these Green functions contain operators defined on up to four sites. Further simplifications can be made if we assume 
that the dominant contributions arise from terms involving at most two nearest- neighbor sites, x and x' . In this 
case, the superscript a in Eq. Q corresponds to a particular neighboring site instead of an average over all the 
nearest-neighbor sites. With this assumption, we obtain 

Dmo{uj) = j^FT{e{t - t'){5J{t),Sj\t')})i, 

Dm,{Lo) = ^^FT{e{t-t'){5J{t),6J'\t')})i, (12) 

where the factor of l/2d arises from the coordination number, a(k) — ^ cos(fci), and 5J and 5 J' are centered on 
two nearest-neighbor sites, x and x', respectively. It is straightforward now to express the Green function in terms of 
these quantities, 

^(k, u:) = ^ \ ,^ ^^^^ ^ (13) 



where the on-site Green function 5*0 is given by 



1, . 1 / (tj + /^)/i +te-i:)TOo(w) -te-|-L»TOo(w) \ I , 

^"^'-"^ 1 -te + L>mo(w) [uj + yi)h +ie - DmQ{u:) ^ ' ^ ' 



and the nearest-neighbor contribution is 



/ 1 + [p - r^Dm^{uj)\l^-' 1 - [p ~ f-ii^mi(c.)]/ri/2-i \ 
^^""^-yx-iv- t-^Dm^{u,)]I^'i:^' l + [p- i-'Dm,{Lu)]I^^ ) ' ^^^^ 



B. The two-site problem: Level operators and resolvents. 

At this stage, solving our problem entails a calculation of the functions DitiqIlu), Dmi{uj) and the parameter p. 
To this end, we express these quantities'^ in terms of correlation functions for the level operators associated with a 
two-site problem. Let us introduce first the single site level operators B{i), F^(i) and D{i) which annihilate empty, 
singly occupied (with spin a) and doubly occupied states, respectively. In terms of these operators, the original 
Hubbard operators can be written as £,a — Fcr and rj^ — aF^^D. As the system can be at a given time in one of 
the possible four states, the level operators satisfy the condition 

Q{{) = i?t(^)i?(^) + J2 Fl{i)F,{i) + D\€)D{€) = 1. (16) 

a 

This restriction can be introduced by adding a Lagrange multiplier term of the form, £oX]t*3(*): to the original 
Hamiltonian. The level operators for the two-site states, are obtained by taking all the combinations of single-site 
level operators. It is convenient to use the following symmetric or antisymmetric combinations corresponding to eight 
fermionic-type states, j 

FB'^s 
FDl 



FD^A 

and eight bosonic-type states. 



-^iF,ix)Bix') + B{x)F^{x')), 
-^{F„{x)D{x')+D{x)F^{x')), 
-^{F,{x)B{x')^B{x)F,{x')), 
-^{F,{x)D{x')-D{x)F,{x')), 



(17) 



BB = B(x)B{x'), DD = D{x)D(x'), 
FF'' = F^{x)F„{x'), 
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FFs ^ ^{F^{x)Fax')+F^{x)FT{x')), 

FFa ^ ^{F^{x)Fi{x')-Fi{x)F^{x')), 

DBs = -^{D{x)B{x')+B{x)D{x')), 
v2 

DBA = -^{D{x)B{x')~B{x)D(x')), (18) 

where x and x' are the positions of the two-sites and the spin index, cr, takes two values, (+1 =t) and (—1 =J,). Note 
that FF'' and FFs correspond to the spin-triplet states, while FFa corresponds to the spin-singlet state. 

We are interested in solving a two-site problem for a cluster embedded in a reservoir, constituted by the rest of the 
system. Formally, the total Hamiltonian can be divided into three parts^ describing the two-site sub-system. Ha, the 
reservoir, TJ^, and their interaction, i?ofl: 

H^Ho + Hn. + HoR. (19) 
To describe the properties of the two-site system, we introduce the resolvent 

^„TrR [(O| $„(t)<i>„ (t0t|o)e-/3gH] 



where the trace is over the degrees of freedom of the reservoir, |0) denotes the vacuum for the two-site problem, and, 
as usual, (3 = l/ksT. Note that $„|0) = and , consequently, i/o_R|0) — 0. The Fourier transform of the resolvent 
can be expressed using the spectral function, cr„„i(cij) = — ilmi?„m, 

Rr,M= I dx ^"'"(^) (21) 

J UJ — X + 10 

We also introduce the auxiliary function CT„m(ijj) = e^^'^anmi^)- Once we know the resolvents, we can express any 
average of operators of the type <i>„m = ^ti^n as 

(*nm) = ^ J doj an„i{uj) (22) 

with Z = J2n I CT„„(w). 

A formal solution of the rcsolvents^^ can be obtained using the equation of motion method. We can write 

RnM = ( , ) , (23) 

where the energy matrix E is determined by the levels of an isolated two-site system and the self-energy E is a measure 
of the effects of the reservoir. Explicitly, we have 

Enrn = (0|z ^ $„ (t)ci>t ^ (i) | Q)^ (24) 

and 

S,„„(cc;) - F.T.{0\e{t - t')5J^^{t)6Jljt')\0)Ri, (25) 

where 'R' indicates that the trace over the reservoir degrees of freedom has been taken, T indicates the irreducible 
part, and (5J$„(t) — i^$„(i) — J^m^nm^mit)- To proceed, we determine the equations of motion for the level 
operators. It is convenient to write first the equations for the single-site operators. 



d_ 

dt 
d 



if^F, = {eo~fi)F,^iBc'^-iacl"^D, (26) 



'di^ = {eo-2fi + U)D-iJ2aF_,c'^, 
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where the arbitrary reference energy Sq will be set to £o = —fJ.- In fact, the lattice and the two-site cluster had differing 
chemical potentials. However, equilibrium between the two-site system and the lattice requires that both have the 
same chemical potential. The equations of motion for the two-site level operators can be determined directly using 
the equations of motion, Eq. I|26|) and Eqs. (|17|l and 1181) . Explicitly, these equations are given in Appendix A. From 
the equations of motion we can extract the energies Enm for the resolvents. Selecting the terms that do not contain 
c"-type operators, that is, the terms that do not depend on the degrees of freedom of the reservoir, we obtain: 

i i 

Edd = 2(60 - 2 fi + U) (27) 



EpBs 


= 2eo- 


Efba 


= 2eo~ 


Ebb 


= 2eo, 


EfFs 


= Ep pa- 


EdBs 


= EdBa 



t 

2eo — 2fi + U, Edbs ffa — Effa dbs — ^ • 

Employing the standard non-crossing approximation^ and using the time derivatives of the $„ operators from Ap- 
pendix A, we compute the self-energies S„,„, Eq. (|25() . of the resolvents. The expressions for the self-energies are 
given in Appendix B. 

It is known^ that inter-site spin fluctuations, which are ignored in the non-crossing approximation, are in fact 
important and cannot be neglected at energy scales on the order of /U . Since our attempt is to develop a top-down 
approach in which on-site physics as well as nearest-neighbour correlations are described accurately, we must include 
spin fluctuations. To overcome this problem, we include the effects of spin fluctuations as a higher order correction 
to the self-energies of the resolvents. Physically, we can understand this correction by observing that the solution for 
singlet and triplet states, FFa and FFs, is sharply peaked at energies separated by 

/oo 
Uj{(TFFs - O-FFa) d^- (28) 
-oo 

In the strong coupling limit, J is of order t'^/U and, consequently, singlet-triplet mixing cannot be ignored. We 
consider J to be the coupling constant of an effective antiferromagnetic interaction which is responsible for the spin 
fluctuations. The corrections to the self-energies given by spin fluctuations are given in Appendix C. The energies 
from Eq. H27|) , together with the self-energies given in Appendix B and the corrections from Appendix C constitute 
the equations necessary for the evaluation of the two-site resolvents. 



C. Self-consistent procedure 



The goal of introducing the two-site resolvents is to express the dynamical corrections Drrio and Drrii, as well as 
the parameter p in terms of quantities associated with the two-site problem as depicted in Fig. To this end, let 
us introduce the quantity Z^^ = Z(<I>}j$„) by 



duj (7$^ (w). 



(29) 



The two-site occupation numbers for each of the 16 states can be written directly in terms of Zg,^ . For example, for 
the singlet and triplet states, we have uffa = Zffa/^ uffs = SZffs/^i respectively. Single-site occupation 
numbers can also be expressed in terms of using Eq. H16() and writing the identity operator on a neighboring site 
as 1' = Q{x'). For example. 



ZJtD = DtDl' = {DB')'^DB' + ^{DF'^)"^ DF'^ + {DD')''DD\ 



(30) 



and as a result 



no 



{D''D) = z-^ 



ZfDs + ^FDa + n^DBs 



'DBa 



Z 



DD 



Similarly, we obtain for the average single occupancy, '^^{FlFa) , 

Uf — Z^^ [3ZfFs + ZfFa + ZfBs + ^FBa + ZfDs 



'FDa\ 



(31) 



(32) 
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FIG. 2: Self-consistent scheme for the computation of the self-energy. 



The electron filling is given by the sum of single occupancy and (twice) the double occupancy, 

n = np ^ 2nD. (33) 

As all the terms in p can be expressed in terms of products of two-site operators, it is straightforward to write this 
parameter as 



3 111 

ZpDs + ZfDa + ^DD + -^ZpFs ^ -^ZfFa ^ -^ZdBs + ^ZdBa 



(34) 



This represents the self-consistency condition for the parameter p and obviates the need to employ a decoupling 
scheme required in the static approximation^fl*^. 

The next step is to determine the dynamical corrections. To this end, we express^ the 'currents' SJ^ in terms of 
two-site level operators. From Eq. lO and Xi ~^ x and a;" —> x' , we find that 

SJa = -iy^anrn^l^rn, 

m,n 

K = -iY.a'^ruK'^rn, (35) 

where $„ represents the complete set of two-sites level operators and the coefficients a„m are given in Appendix D. 
As dJ' is obtained from 6J by exchanging the positions x and x' , the coefficients a'j„j will be identical with anm up 
to a sign that depends on the symmetry properties of the states '$m' and '<&„' under the exchange of x and x' . Let 
us denote by £ the operator that produces the exchange, 

£0{x,x') ^ 0{x',x), (36) 

where 0{x, x') is an arbitrary operator defined on the two-site cluster. The symmetry properties of the level operators 
are given by 

EFBs = FBs, £FBa = -FBa, SFDs = FDs, £FDa^-FDa, 

EFFs = -FFs, EFF" = -FF" , EFFa=FFa, (37) 

EDBs = DBs, EDBa = -DBa, EBB = BB, EDD = DD. 

Consequently, = a„„i if $„ and $,„ have the same symmetry and a'^.^ = —anm otherwise. The self-energy 
contributions can be calculated only approximately. We will use the non-crossing approximation^-'^ which has been 
proven to be effective in solving problems with local correlations. Defining <f>mn = ^n^mi we have 

Trn mrnnml,Ai'y-^"n (38) 
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~ TrR{Q\<^^ml,{t')\Q)TrR{Q\^^,{t')<^i{t + z/?)|0). 

Consequently, we obtain 

($m„(0*L'„'(i')> ~ I /^^ e-''^(*-*') j dx aw(w + x)a,,,n{x), (39) 

which is the formulation of the non-crossing approximation that we will use systematically in our calculations. Intro- 
ducing the expansion Eq. (|35|) in the expression for the dynamical corrections, Eq. p2|l . and using the non-crossing 
approximation formula( Eq. H39() ). we obtain 

r-i / \ 1 fill \ ^ it ^rnm' i-^^^n' ni-^ ) ^ mm' i'^') ^ n' ni*^ ) r Ar^\ 

Dmoiu;) = — dxdx' a™a*,,„ — (40) 



n.m.n' ,m' 



2dZ J ^ , wx + x' + i5 

and 



DmA-) = ^ dxdx Y: .^x + x' + ^S ^^^^ 



n,m,n' ,m' 



Detailed expressions for Dtuq and Dmi are given in Appendix D. 

This concludes the process of writing the components of the Green function Eq. jSjl that cannot be expressed in 
terms of the Hubbard basis, in terms of spectral functions for the two-site system. However, for a fully self-consistent 
calculation, we need to determine the reservoir spectral functions p±{uj) that enter into the formulas for the self- 
energies of the resolvents (see Appendix B). The reservoir consists of the full system from which the two sites x and x' 
have been excluded. Our goal is to determine the Green functions G = ((c", cj.")) and G" = ((c", cj." )). We introduce 
the two-site full propagator [5] = ((5', ^f^)) with — {tp^ {x)^p\x')) and the irreducible propagator 

[S] (42) 



^S' s ^ 

where S = ((-(/'", -0^")) and 5" = )). From standard scattering theory, we have that 

[5] = [5o] + [5o][F][5][^][5] (43) 

where 
and 

The solution of these equations is 

s = v'^So^ - {s- s's-^s')-^)v~\ 

S' = V-^ {tV + S-^S'{S - S'S-^S')-^)V-^. (46) 
The irreducible Green functions used in the evaluation of the self-energies of the resolvents will be given by 

G = Sii + Si2 + S21 + S'22, 

G' = S[,+S[^ + S'^, + S'^„ (47) 
and the corresponding spectral functions are 

p+iio) - --Im(G(^) + G'(a;)), 

IT 

p-{uj) = --Im(G(w)-G'(w)). (48) 
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III. RESULTS 



We now have all the ingredients necessary for the implementation of the self-consistent procedure shown in Fig. jSJ- 
Starting with an initial guess for the spectral functions p± describing the properties of the environment, we solve the 
two-site problem iteratively using the expression Eq. H23|) for the resolvents, the energies Eq. (|27|l and the self-energies 
from Appendix B with the spin fluctuation corrections from Appendix C. While in principle the cluster and the lattice 
can have different chemical potentials as in the previous work on the two-site cluster—, we have used the more physical 
restraint that both the cluster and the lattice must have the same chemical potential. However, as a result of Eq. 
(|33|l . the filling in the cluster and the lattice will be different. By symmetry, however, at half-filling both the cluster 
and the lattice will have the same filling. All the frequency-dependent functions are discretized on a grid of TV = 8192 
points from LUmin — —20t to LOmax — 20i. To increase the computational speed, we performed all the convolutions 
involved in the calculation of the self-energies using a fast Fourier transform algorithm. The procedure converges for 
temperatures above T = 0.02t at finite doping and T = O.OSt at half-filling, although convergence problems occurred 
below T = O.lt. Once we computed the resolvents, the mean- field parameter p can be determined using Eq. H34() . as 
well as the dynamical corrections DrriQ and Dmi from Appendix D. The mean-field parameter e = ~ (^"^^) 

can be expressed in terms of the Green function using the general self-consistency condition. 



i2ny 

X (^~~^ ImS'm„(k, w). 

(49) 

Within the grand-canonical ensemble, the chemical potential, fj,, is determined by the self-consistent solution to 

n = 2 {((^^ + 2^Tj^ + r]r]^)) (50) 

We imposed the constraint that the chemical from Eq. H5U|) also equal that for the cluster. We then determined the 
full Green function S'(k, w) using Eqs. (|13I15|) . Next, new spectral functions p± are determined using Eq. H48|) and 
the whole procedure (see Fig. |(5J is repeated until full convergence is reached. 



A. Spectral Function at Half-filling: Mott Insulator 

Before we analyze the doped case, we first review2S the properties of the charge vacuum that determines the 
insulating behavior at half-filling. To reiterate, there are two distinct routes to the insulating state at half-filling. 
For U ^ t, the charge and spin degrees of freedom are decoupled and the system is an insulator for temperatures 
smaller than Tq ~ U. The spins are coupled due to the super-exchange interaction, |J| At^ /U. It is this spin 
exchange interaction that gives rise to local antiferromagnetic fluctuations and eventually ordering at T = 0. This 
is the antiferromagnetic Mott insulating state. In general, in the weak coupling regime, a metal-insulator transition 
occurs as a consequence of the Brillouin zone folding generated by magnetic or charge ordering and the corresponding 
gap is essentially related to antiferromagnetism or some type of charge density wave. This type of transition is referred 
to as a Slater transition and the corresponding insulating state should not be confused with the Mott insulator. Such 
a regime can be successfully described by conventional many-body approaches. 

Shown in Fig. Q is the total electron spectral function, 

v4(k, uj) = --Im {5n(k, uj) + 25i2(k, uj) + S22{K ^)} (51) 

TT 

for the 2D half-filled Hubbard model with U — 8t and T = 0.15t. Clearly visible are the upper and lower Hubbard 
bands with an energy gap of order U and flatness of the band near the (tt, 0) point. The chemical potential {lo — 0) 
lies in the middle of the gap and hence the system is an insulator. However, no symmetry is broken as is evident 
because the periodicity is 27r rather than tt as would be the case if the Brillouin zone had doubled. In addition, spin 
and charge are not fractionalized^^. The insulating behavior arises because the charge gap has splintered the spectral 
weight of each k-state into 'bonding' and 'antibonding' pieces. Consequently, there is a fundamental breakdown of 
what is meant by an electronic state. In fact, the electronic states themselves have fractionalized. To make contact 
with the real-space picture shown in Fig. J^l, we note that the PES and IPES spectra are determined by the spectral 
weight in the lower and upper Hubbard bands, respectively. In general, the upper and lower Hubbard bands carry 
total spectral weight n/2 and 1 — ri/2, respectively which of course reduces to 1/2 at half-filling. However, A(k,uj) is 
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FIG. 3: Momentum and energy dependence of the electron spectral function for a half-filled 2D system with U = 8t and 
T = 0.15t. From top to bottom, the curves correspond to {k^, ky) = (0, 0) (tt, tt) (tt, 0) — > (0, 0). 




FIG. 4; Occupancy of each momentum state with kx — ky for the spectral function shown in Fig.Q. The occupancy is peaked 
at (0, 0) and decreases continuously to a minimum value at (tt, tt) without the discontinuity indicative of a Fermi liquid. 



strongly momentum dependent as illustrated in Fig. (PJ. For the lower Hubbard band, the maximum in the spectral 
weight is peaked at (0, 0) and decreases as the (tt, 0) point is reached and becomes vanishingly small at (tt, tt). In fact, 
the states at (tt, tt) in the lower Hubbard band carry almost none of the spectral weight. However, the decrease in the 
occupancy of each k-state is a continuous function, as depicted in Fig. |0J, rather than a discontinuous one as would 
be the case in a Fermi liquid with a well-defined Fermi surface. In the strict sense, the discontinuity in rik in a Fermi 
liquid occurs at T = 0. Although we cannot reach T = in our approach, we find no indication that a discontinuity 
develops in n|„ecfc as the temperature decreases. In fact, the continuous behavior we have obtained here is consistent 
with the exact resul'b^ 

"k = ^ + 2^(S,.S,+,-i) (52) 

for the occupancy in each k— state projected into the LHB of a half-filled Hubbard model. Note all the k— dependence 
is determined by the single particle energies, 

ek = -i^e''"(^'-^*+^) (53) 

s 

For a paramagnetic state, the term in the angle brackets is exactly, —1/4. Hence, is continuously decreasing 
function from (0, 0) to (tt, tt) as found here for U = 8t. 
Although rik < 1, each state satisfies the sum rule. 



/oo 
A{k,uj)duj = 1. 
-oo 



(54) 
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In a band insulator the upper cutoff on the energy is simply the chemical potential. This does not mean that in a 
band insulator or in a Fermi liquid, a broad incoherent background cannot be present which extends to high energies. 
In fact, in a Fermi liquid, can be less than unity. However, in Fermi liquids, a coherent quasi-particle peak always 
exists regardless of the momentum. Hence, the criterion for occupancy of a single particle state is simply whether or 
not the coherent peak lies above or below the chemical potential. In a Mott insulator not only is the spectral weight 
split over an energy scale of U but there is an absence of coherent quasi-particles as evidenced by the broad spectral 
features. Hence, there is no sharp criterion for unit occupancy of a single particle state. The broadness of the spectral 
features stems from the local correlations on neighbouring sites not the Mott gap itself. Without the dynamical 
corrections, the spectral function would simply be a sum of (5— function peaks at the lower and upper Hubbard bands. 

The bifurcation of the spectral weight of each k— state above and below the charge gap can be modeled as follows. 
Consider for example 7^^ = ""kCkcr + ^k?/ko-' with coefficients and Wk are determined by the projection of the 
spectral function onto the lower and upper Hubbard bands, respectively. Hence, the antisymmetrized state formed 
from such excitations 

|MI)=^(-lf n 7kT7ljO) (55) 

P keFBZ 

is a candidate for describing the elusive paramagnetic Mott insulator. Provided magnetic frustration is present so 
that ordering is pre-empted, Eq. 1(531) should be the T = Mott insulating state. The sum over all permutations, P, 
is necessary as the 7kCT operators obey the non-fermionic commutation relations, 

Hv^ll} = E^^^^'^'^^^'K^^ - ^'^^^K4l (56) 

i 

and {7kCT) Iqa} — {ita' 7qo"} — 0. When Uk — Vk, lt\lki\^) generates the completely doubly occupied state. However, 
because the right-hand side of the anticommutation relation in Eq. H5t)|l is identically zero for Uk = v^, the zero state 
results upon summation over all permutations. Consequently, Eq. H55|l completely projects out the fully doubly 
occupied state. 

1. Singlet Formation: Local Antiferromagnetism 

A crucial test of the correctness of the method we have used here is whether or not short-range antiferromagnetic 
correlations are present at low temperatures. Such correlations do not signal that long-range magnetic order obtains 
at T = but rather that the ground state at T = is a liquid of nearest-neighbour singlet states as in the resonating 
valence bond (RVB) state proposed by AndersonSi. We are able with our two-site formalism to probe the existence 
of local magnetic order by computing the nearest-neighbour singlet and triplet occupation numbers, uffa — Zffa/^ 
and uffs — '^Zffs/Z, respectively. From Fig. (O we find that, at high temperatures, triplet excitations dominate. 
However, this trend is reversed below some temperature and the singlet occupancy becomes of order unity. Hence, 
the low-temperature properties of the insulating state we have computed here are consistent with a liquid of nearest- 
neighbour singlet states as in the RVB state. In fact, at T = the liquid state we have found here persists because 
we have imposed a paramagnetic solution. 

A final diagnostic of the insulating state we have found here is the behavior of the effective exchange interaction as 
a function of U. In the Mott state, a super-exchange interaction is self-generated which should scale as 1/U. It is this 
exchange interaction that sets the scale for the Neel temperature. Using Eq. (|28|) . we computed the effective exchange 
interaction shown in Fig. © for both ID and 2D. Note first that J is always positive as a consequence of the fact 
that the singlet state is lower in energy than the triplet. This is a further indication of antiferromagnetic correlations 
in the ground state. As expected, J is well approximated by /U in the strong-coupling regime. However, as U 
decreases, deviations from this behavior are observed. 

B. Doped Mott Insulators 

1. Chemical Potential 

Two scenarios are possible for the doping dependence of the chemical potential: 1) the chemical potential remains 
pinned and mid-gap states are generated by some physical mechanism, or 2) the chemical potential jumps to the top 
of the lower Hubbard band (LHB) or the bottom of the upper Hubbard band (UHB) upon hole or electron doping, 
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FIG. 5: Singlet (wff^) and triplet {npFs) occupation numbers as a function of temperature for U = 8t in two dimensions. 
The fact that tiff a > ^ffs as T — > is consistent with local antiferromagnetic order. 
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FIG. 6: Effective exchange interaction coupling constant J as a function of U/t for T — 0.15t. 



respectively. Our results shown in Fig. (O demonstrate that the chemical potential jumps upon hole or electron 
doping, indicating an absence of mid-gap states. The Magnitude of the jump is set by the Mott gap which is fully 
developed at T = 0. While at some finite temperature the chemical potential may appear to evolve smoothly, at 
T = 0, the chemical potential should jump discontinuously. That is, we find that doping a Mott insulator leads to 
a continuous depletion of the spectral weight in the first Brillouin zone. Our finding is consistent with the exact 
result for the ID Hubbard modelS as well as Quantum Monte-Carlo simulationgiSii^ in 2D. However, in d = oo, the 
chemical potential exhibits^^ a jump but one that is smaller than the gap. Hence, mid-gap states are generated^^. 
It would appear then that d = oo is vastly different from an actual finite-dimensional system and the 2D Hubbard 
model is quite similar to its ID counterpart at least as the chemical potential is concerned. The possible source of the 
discrepancy is the form of the self energy. A chemical potential jump requires a large imaginary part of the self-energy 
at the chemical potential, thereby indicating an absence of well-defined quasi-particles. Mid-gap states are resonance 
states and hence are reminiscent of the Brinkman/Rice^ mechanism for the insulator-metal transition in the doped 
Mott state. The inset in Fig. ^ indicates that in the underdoped regime is large and non-zero at the Fermi 
energy. Such behavior points to an absence of well-defined quasiparticles. In the overdoped regime, the characteristic 
dependence appears, indicative of a Fermi liquid. Consequently, the method we use here is capable of recovering 
Fermi liquid theory in the overdoped regime. 

Experimentally, whether the chemical potential is pinned or moves upon doping appears to be cuprate dependent. 
For example, in La2-xSra;Cu04^* (LSCO), the chemical potential remains pinned roughly at OAeV above the top of the 
LHB, while for Nd2-^Ce:,Cu04 (NDCO)M, BiaSraCai-^R^CuaOg+y (BSC0)SMiiS2iS, and Na-doped CazCuOzCla^ 
(CACLO) the chemical potential jumps upon doping by an amount in accordance with half the Mott gap and scales 
roughly as 6^ as obtained here. Because stripes or macroscopic phase separation require the chemical potential to be 
pinned, they have been invoked?^ to explain the origin of mid-gap states in LSCO. The pseudogap in the underdoped 
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FIG. 7: Doping dependence of the chemical potential in the 2D Hubbard model computed using the local cluster approach for 
T = 0.15t (dashed line) and T = 0.07t (solid line). The inset shows the imaginary part of the self energy evaluated at a Fermi 
momentum (0.3, 2.10) for n = 0.97, (0.3, 1.84) for n = 0.8 and (0.3, 1.06) for n = 0.3. 

cuprates has also been attributed^ to stripes. However, because A/i ^ for most of the cuprates, for example, 
NDCO, BSCO, and CACLO, if the pseudogap has a universal origin, stripes are not its cause. The precise origin of 
the pseudogap will be discussed extensively in a later section. 

2. Spectral Function 

Shown in Fig. ||SJ| is the electron spectral function at high temperature, T = .25t for n = 0.97, n = 0.90, n = 0.80, 
and n = 0.60. Several features are evident: 1) the chemical potential moves further into the LHB as the filling 
decreases, 2) no coherent peaks exist near the chemical potential in the lightly doped regime, .9 < n < 1, 3) in the 
dense or weakly interacting regime, sharper features appear, 4) each state in the FBZ has spectral weight both above 
and below the chemical potential as dictated by Mottness, 5) the Mott gap remains intact but moves to higher energy 
as the doping increases, and 6) at {tt,tt), the UHB carries most of the spectral weight regardless of the filling. In the 
underdoped regime, the characteristic width of each k state is of order t and even much larger near (tt, 0). Such broad 
spectral features in the underdoped regime are seen experimentall}i^ and arise in this context because SJE(eF) 9^ 
as shown in Fig. 0. As a consequence, there is no sharp criterion for unit occupancy of each state in the FBZ. 
Because the total spectral weight of each k state is unity, however, and each state lives both below and above the 
chemical potential, the charge carried by the piece of the state lying below the chemical potential is less than unity. 
That is, each electronic state is fractionalized. In the heavily overdoped regime, the splitting of the spectral weight 
above and below the Mott gap, is highly suppressed. As illustrated in Fig. ((SJ, most of the spectral weight resides in 
the LHB for a filling of n = 0.60. As a consequence, Mottness vanishes in the overdoped regime. In addition, in the 
heavily overdoped regime, acquires the characteristic dependence indicative of a Fermi liquid. Hence, with our 
method we are able to recover the key characteristics to the transition to the traditional Fermi liquid state, namely 
1) a vanishing of spectral weight in the UHB and 2) « near the chemical potential. 

Does new physics emerge at low temperatures? Fig. depicts the spectral function computed at T = 0.07t. 
At this relatively low temperature, two new features emerge. First, in the underdoped regime, the spectral weight 
appears to be suppressed at the chemical potential. Whether this gives rise to a pseudogap will be investigated in the 
next section. Second, at n « 0.8 especially in the vicinity of the (tt, 0) point, the band becomes almost dispersionless 
and seems to split into two sub-bands. Instead of a strong coherence peak, at the chemical potential we observe a 
weak maximum adjacent to a region with depleted spectral weight which forms a 'channel' immediately above the 
chemical potential. To understand the importance of the low-temperature features, it is expedient to compute the 
density of states. 
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FIG. 9: Doping dependence of the spectral function in the 2D Hubbard model computed using the local cluster approach for 
T = O.On and fillings of n = 0.97, n = 0.90, 0.80, and n = 0.60. 
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FIG. 10: Density of single particle states for T = 0.25t and T = Q.Qlt, U = 8t for the fillings shown. No pseudogap exists 
at high temperature. At low T and low doping levels, a pseudogap emerges at the chemical potential but moves above it at 
an intermediate doping level. In the overdoped regime, the pseudogap vanishes entirely and a weakly-interacting system is 
recovered. 



3. Pseudogap without preformed pairs or global symmetry breaking 

To investigate the possibility of a pseudogap^^ in the lightly doped regime, we integrate the spectral function over 
momentum to obtain the single particle density of states (DOS) at high and low temperatures. Displayed in Fig. 
(|10|l is the DOS for T = 0.25t and T — 0.07i for several fillings. As is evident, no local minimum of DOS exists 
at the chemical potential at high temperature, T = 0.25t. Features which emerge even at high temperature are the 
reshuffling of spectral weight from above the charge gap to below as the filling is changed and also a movement of 
the Mott gap to higher energies. Note that even at n = 0.30 the Mott gap is still present, though almost all of 
the spectral weight now resides in the LHB which closely resembles the non-interacting density of states. This is 
further evidence that we correctly recover Fermi liquid theory as n — > 0. What about low temperature? The lower 
panel of Fig. (|10|l demonstrates that a pseudogap forms in the DOS for 5 « 0. The vertical line at indicates 
that the pseudogap occurs precisely at the chemical potential. Similar qualitative results based on a cluster method 
have been obtained by Maier, et. aliS&, except their pseudogap is slightly displaced above Ep. In contrast, in the 
analysis of Haule, et. al&, the DOS has a negative slope through Ep (as dictated by the proximity to the Mott gap) 
but never acquires a local minimum at Ep indicative of a true pseudogap. Because the pseudogap exists below some 
characteristic temperature and vanishes at higher doping, the result obtained here is non-trivial and highly reminiscent 
of experimental observations in the cuprates^. What is its origin? Incoherence (3E ^ at Ep) is central, though it is 
not a sufficient condition-7{ see Fig. 2) for a pseudogap. From Fig. l(TU|l . the pseudogap remains intact up to (5 = 0.20 
but simply moves to higher energies as does the Mott gap. This is telling because in d = oo, no^ pseudogap exists 
but a Mott gap is present. Absent from d = oo but present in any lattice of finite connectivity are true short-range 
correlations. We argue then that the pseudogap is the nearest-neighbour analogue of the on-site generated Mott gap. 
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FIG. 11: Density of single particle states for T — O.lt, for n — 0.95 and three values of the on-site interaction: U = 12t, U = 8t, 
and U = 4t. The gradual vanishing of the pseudogap as U increases offers direct confirmation that the energy scale for the 
pseudogap is set by /U. 



The energy scale for nearest-neighbour interactions scales as /U . Hence, if our hypothesis is correct, we expect 
the pseudogap to diminish as U increases. The evolution shovifn in Fig. Hll|l indeed demonstrates that, at finite 
temperature, the pseudogap does vanish as U increases. Hence, we can assert with certainty that correlations on 
neighbouring sites do in fact create a depletion in the density of states. The energy scale /U is typically associated 
with antiferromagnetic spin fluctuations. To explore whether such processes have the right doping dependence to 
explain the origin of the pseudogap, we display in Fig. (|13|l the x— dependence of J computed from Eq. H28|l. As is 
evident, J is only weakly doping dependent in the underdoped regime and hence lacks the strong doping dependence 
needed to explain the pseudogap. This trend is consistent with that of Jarrell and co-workers^ who have observed 
that the pseudogap persists even if antiferromagnetism is killed. Fig. (|13|) also indicates that Jcg computed as the 
energy difference between the nearest-neighbour singlet and triplet states vanishes at a; = .8. This is not an accident. 
Nearest-neighbour spin fluctuaions should desist when no nearest-neighbour sites remain singly occupied. On average, 
this obtains at a filling of 1/5 or a; = .8, precisely the doping value found here. Hence, it is not a surprise that the 
doping dependence of J is weak in the underdoped regime. 

What then is the cause of the physics of the pseudogap? Any two-step process involving the UHB scales as /U . 
Consider the explicit three-site terms that appear away from half-filling 



in an expansion in t/U , where the tilde represents full projection onto the LHB. It is precisely these terms that thwart 
the equivalence between the so-called t — J and Hubbard models in the large U limit away from half-filling^^. This 
term represents the motion (strictly in the LHB) of a hole in a spin background. We argue that such terms are 
involved in the pseudogap. The mechanism is as follows. Consider placing a single hole in a Mott insulator. Unlike 
a site neighbouring the hole, a singly-occupied site two lattice sites away must temporarily doubly occupy one of its 
neighbours if it is to move to the hole. For this to be possible, the electrons on neighbouring sites must have opposite 
spins. The matrix element for such a two-step process is /U and described by the three-site terms written above. For 
sites with the incorrect spin alignment, a local spin fluctuation must obtain for the three-site hopping to occur. The 
energy barrier for this process is /U . It is from those local three-site configurations in which the spins are incorrectly 
aligned that the pseudogap arises as illustrated in Fig. (|12ll . Simply invoking spin fluctuations is insufficient to explain 
the origin of the pseudogap as spin fluctuations alone cannot give rise to transport. However, spin fluctuations can 
make it impossible for an electron two sites away from a hole to transport. Hence, spin fluctuations in the context of 
three-site hopping can overcome the local spin blockade (or spin gap) that exists in doped Mott insulators. As this 
effect is entirely local, the pseudogap is the nearest-neighbour analogue of the Mott gap: neighbouring sites with a 
parallel arrangement of the spins experience an energy barrier equal to /U for charge transport. Can the doping 
value at which this process vanishes be estimated? On this account, the pseudogap should be related to the joint 
probability that a neighbouring three-site configuration consiting of a hole and two sites with spin parallel electrons 
exists. The minimum constraint however is simply that each site has on average one hole as its immediate neighbour, 




(57) 
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FIG. 12: Local three-site configurations in which spin-blocking leads to a pseudogap. In the upper state, transport directly 
between the up spin two sites away from the hole and the hole is not possible. Only hole transport via a two-step process is 
possible. As the amplitude for hole transport is a superposition of all such processes, a pseudogap develops. The only way to 
overcome the spin-blockade is for the spin neighbouring the hole to flip. This process costs an energy /U. Once the spin is 
flipped, there is no barrier for transport. 
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FIG. 13; Effective nearest-neighbour J as a function of filling. J vanishes when four out of five sites is empty. Beyond this 
concentration, nearest-neighbour spin fiuctuations are not possible. 



roughly x ~ 0.25 for a square lattice. Hence, the pseudogap is of the form /U P{x), where P{x) determines the 
probability that hole transport involves double occupancy and consequently, is a steadily decreasing function of x 
vanishing at Xcrit- From the estimate given above, it is likely that Xcrit is closer to the end of the superconducting 
dome than it is to optimal doping. 

It is commort3S-Ifi, to approximate the three-site terms considered here as xKJ, where K is treated as the lattice 
connectivity. Consequently, the effective nearest-neighbour exchange interaction is doping dependent and given by 
Jc« — — xK) which vanishes at Xait — 1/ K k, .25^^ for a square lattice. However, as Fig. demonstrates, Jeff 
vanishes at a; = .8 not x = .25. It is likely that the discrepancy found here arises from the fact that the 3-body terms 
in Eq. (|57|l cannot, in any real sense, be rewritten as an effective spin-exchange. In fact, &t x — .8, every site has on 
average four neighbouring holes. At this concentration, nearest-neighbour spin fluctuations are not possible; hence, 
Jeff should vanish. 
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4-. Heat Capacity 

Evidence for the pseudogap is also found from the heat capacitj*£^. We computed the heat capacity numericaUy 
from the internal energy 

,5S, 

The energy per site is the sum of the kinetic term and the interaction term 

^^-2t{clc^a) + U{n^^n,i), (59) 
where t = 2dto {t = Ato in two dimensions). The double occupancy can be expressed as 

{n,^n,g) ^ (4^r],„) = | - iVi'jcja)- (60) 
The correlations can be expressed in terms of the Green function as 

En f rill) f rl^k —1 —1 

= oU+ — —^(1 - f{Lo))[2ta{k)—Im{Sn + 2Su + S22) ~ U—Im{S^2 + ^22)]. (61) 

Equivalently, we can use the equation of motion for the Green function ((?yi(''"), cj(T'))) and we obtain an alternate 
expression 

J ^{^-fH){u;-ta{k))^ImiS^, +23,2 + 822) (62) 

for the energy per particle. We found that the difference between Eqs. (|61|l and 162() wqs within our numerical 
errors. Gonsequently, in our final calculations of the heat capacity shown in Fig. (|14|) . we averagedthe two results. 
In the ID and 2D Hubbard models at half-filling22i, two peaks exist in the heat capacity. The high temperature peak 
corresponds to charge excitations and the low-tempeature peak to spin physics. As is evident, the same separation of 
energy scales persists even in the doped case. However, for n < 0.9, we find that the spin peak vanishes and merges by 
the charge excitation spectrum. This dramatic change represents a possible termination of Mott-dominated physics 
and the onset of more Fermi liquid behavior. Bonca and PrelovsekSi observed the identical trend in their exact 
diagonalisation study of of a 4 x 4 system. This agreement lends further credence to our method. Another trend 
evident from Fig. (|14|) is that the extrapolated T — value of the heat capacity in the underdopcd regime, n > 0.9 
is lower than that a,t n ~ 0.85. This represents a loss of spectral weight at low energies as would be the case once a 
pseudogap opens. Hence, the thermodynamics also corroborate the existence of a pseudogap. 



5. Possible time-reversal symmetry-breaking 



We have shown that in general, a pseudogap exists in a doped Mott insulator without invoking symmetry breaking of 
any kind. Nonetheless, we entertain the possibility that in analogy with antiferromagnetism and the Mott gap, perhaps 
some broken symmetry state obtains at lower temperature as a result of the pseudogap found here. Our argument 
should be construed as a conjecture and hence is entirely speculative. It by no means underlies the calculations we 
presented here. Loram and colleagues2^ have argued that for x < Xcrit, a- glassy phase with an Edwards- Anderson 
order parameter obtains. While the experimental evidence for a glassy phase extending to Xcrit is not clear, recent 
circular dichroism experiments^^ point to time-reversal as the relevant symmetry that is broken in the pseudogap 
phase. However, this symmetry is broken only along the (7r,0) and (0,7r) directions and not along (tt,tt). Should 
these results endure, they will provide a benchmark for measuring the validity of the numerous proposals for the 
pseudogap£2i2ii22iSi24. In the context of the view put forth here, we must determine how purely nearest-neighbour 
correlations can give rise to a breaking of time-reversal symmetry only along the canonical x and y axes but not along 
X — y. Consider the three-body term in Eq. (|57|l . This term generates correlated motion of a hole among nearest- 
neighbour sites, that is, local currents. In analogy with the local moments that order antiferromagnetically T = as a 
result of the Mott gap, we propose that the currents may order in the pseudogap phase below T*. Experimentally^^, 
translational symmetry is preserved in the pseudogap phase. Hence, staggered orbital currents are not possible as 
they automaticaly result in a doubling of the unit cell within a single-band model^^. Further, experimentally22i2SiSi, 
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FIG. 14: Heat capacity computed for three fillings computed by numerically difi^erentiating the internal energy obtained from 
the average of Eqs. I61II and 16211 . For all three fillings shown, U — 8t. 




FIG. 15: Current pattern for the motion a hole in a doped Mott insulator that preserves t-reversal symmetry along x — ±j/ but 
violates it along the canonical x and y axes. Each lattice represents a copper site. Hence, translational symmetry is preserved. 
The diagonal current is chosen such that the net current in each plaquette vanishes. 



there is no evidence that physics beyond a single-band is relevant to the cuprates. Consequently, any current pattern 
must preserve translational symmetry within a single band model. Only one option remains: the currents order below 
some characteristic temperature, T*, along the canonical x and y axes. To ensure that the net current along x — y 
vanishes, a compensating diagonal current must be present as depicted in Fig. H15|) . This current pattern can be 
obtained from the most recently proposed pattern of Simon and Varmai^ by simply integrating out the oxygen sites. 
Hence, despite claims to the contrary^^, it is entirely possible to generate a translationally invariant current pattern 
within a 1-band model that is consistent with the experimental observations. In the corrected pattern of Varma and 
Simon^, the oxygen sites do nothing except produce a diagonal current which ensures that the total current in each 
plaquette vanishes. That the oxygen sites can be integrated out is certainly consistent with the now well accepted 
work of Zhang and Rice*'^. Nonetheless, our work does not hinge on the current pattern shown in Fig. (|15|l being the 
origin of the pseudogap. However, insofar as such a pattern obtains entirely from local nearest-neighbour physics, it 
is consistent with our finding that the pseudogap is the second harmonic of the Mott gap. Should further experiments 
confirm the presence of t-violation, then a more microscopic investigation of the origin of the current pattern shown 
in Fig. H15|l will be warranted. Of course, a current pattern of the type proposed here can only be obtained (if at all) 
from a Hubbard model if next-nearest neighbour hopping is included. 



22 



0.8 



0.2 











\ 




\ / 


\N / 


w,, \ / 
^ \' 


^ / 

\ A 






' \ 


/ V N 


\ 

^ \ 
^ \ 


/ ^ 

/ 


^ " " \ 


/ 



0.5 1 1.5 2 

Filling n 

FIG. 16: High {Wh) and low (Wl) spectral weight as a function of filling. Wni is the spectral weight in the non-interacting 
system. 



6. Spectral Weight Transfer 



To quantify the spectral weight transfer evident in Fig. (|10|l . we compute the high and low spectral weight by 
integrating the DOS from a value of the energy inside the Mott 'gap' which minimizes the DOS to oo (— oo for 
electron doping) and from fx to that fixed energy, respectively. The results shown in Fig. Hlfifl demonstrate that 
the initial spectral weight in the UHB which is 1/2 at n = 1 all moves to low energies as the filling decreases, as 
in observed experimentally^Siffi. The integrated spectral weight has been normalized per spin. The same is true for 
electron doping (n > 1). Further, the curvature of the low energy spectral weight is positive as a function of doping 
in agreement with earlier resulta^S on the ID Hubbard model. This signifies that the integrated low-energy spectral 
weight increases faster than 2x. The additional low-energy spectral weight above that dictated by state counting 
(see Fig. arises from virtual excitations between the lower and upper Hubbard bands. Such virtual transitions 
arise from the three-site terms discussed previously.. That the number of low energy degrees of freedom arise from 
such high energy processes further attests to the inseparability of the low and high energy degrees of freedom in a 
strongly correlated system. For contrast, the spectral weight transfer for a non-interacting system (Wni) is shown as 
well. This behavior is expected for a doped band insulator. That Mottness leads to such a drastic deviation from the 
non-interacting result is a direct consequence of state fractionalization. That is, each state has spectral weight both 
above and below the chemical potential (see Fig. 



7. Hall Coefficient 



Experimentally, the Hall coefScient in the cuprates is in general positive in the lightly hole-doped regime, scales 
as l/x in the vicinity of half- filling but falls off faster than 1/x for a; « .1 and in some instances changes sign2^iSSi2& 
typically around x = .25. Although the Hall coefficient is in general temperature dependent as emphasized extensively 
by AndersorS, the Takagi, et. al. experiments®''^ indicate that the zero-crossing doping level in La2_a;Srj.Cu04 
(LSCO) is only weakly temperature dependent. Consequently, we will focus solely on the doping dependence of the 
Hall coefficient since the existence of the zero-crossing is only weakly doping dependent. Nonetheless, because the 
sign change is not universally observed in all the cuprates, the general conditions under which a sign change of the 
Hall coefficient should be observed in a doped Mott insulator have not been formulated. In addition, there have 
been numerous theoretical treatments of the Hall coefficient. For example, perturbative schemeSi2^*22i22iSS lead to a 
sign change of Rh and hence offer a possible explanation for the deviation from l/x. However, because perturbation 
theory is constrained by Luttinger's theoremiSi to yield a Fermi surface occupying half the first Brillouin zone (FEZ) 
at half-filling, such approaches fail to recover the experimentally observed^ divergence of Rh at half-filling. In 
strong-coupling calculations, some have obtained a sign changeS^i^SiSi while othersS^ predict that Rh < for all hole 
dopings. In addition, others2Si2I have reached the counterintuitive conclusion that Rh does change sign, but the 
Fermi surface is closed for all a; > 0. In such studies, it was assumed that the doped Mott insulator is described by 
doping the diamond-shaped Fermi surface of the weakly-interacting system, an assumption clearly not borne out by 
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FIG. 17: (Color) Spectral function in the first Brillouin zone evaluated at the chemical potential for a filling of n = 0.97 (top 
panel) and n = 0.30 (bottom panel) with U — 8t and T — .25t. Color scheme: blue= minimum (zero) spectral weight, red= 
maximum, yellow= intermediate values. In the underdoped regime, the spectral function has broad features at the Fermi level. 
A sharp Fermi surface emerges in the overdoped regime. 



experimen1iSi22i2£. 

On simple grounds, however, the general doping dependence of the Hall coefRcient can be easily deduced. Consider 
a lightly-doped Mott insulator in which the Hall coefficient is initially positive. In the heavily overdoped regime 
where the system is weakly interacting, the Fermi surface must be closed and hence electron-like; thus Rh < 0. We 
can deduce the doping level at which the transition from an open to a closed Fermi surface occurs by appealing to 
the spectral function of a Mott insulator. As Figs. (O and I0J illustrate, at half-filling every k— state in the FEZ 
has some spectral weight. Because the chemical potential (see Fig. {T))) simply moves down through the LHB upon 
doping, hole-doping simply depletes the spectral weight in the LHB. When half the spectral weight in the LHB is 
removed, the surface separating occupied from empty states must have zero curvature. As this surface defines the 
Fermi surface, then its curvature should be related to the sign of the Hall coefficient. Consequently, the critical doping 
level at which Rh = is determined by the doping level at which half the spectral weight in the LHB is depleted. 
Consider the U — oo limit in which the spectral function is momentum independent. Each state in the LHB carries 
the weight 1 — n/2. Consequently, the fraction of the spectral weight depleted upon hole doping is n/2/{\ — n/2). 
When this quantity equals 1/2, Rh should vanish. The solution to 

— ^ = — lower limit on the filling at which Rh = Q (63) 

is n = 2/3 which is a strict lower bound for n^ and is precisely what simulations as well as complicated series expansions 
in the infinite-L'' limi-feS^t^^ obtain. Consequently, if a Mott insulator possess a sign change in the Hall coefficient, it 
must occur for x < .333. 

Shown in Fig. (|17l) is the spectral function in the FBZ evaluated at the chemical potential for U = 8t and T = .25t. 
The upper panel corresponds to n = 0.97 and the lower to n = .3. As is clear, in the lightly-doped regime, the Fermi 
surface is hole-like and the spectral features are broad indicating an absence of well-defined quasi-particles indicative 
of an incoherent metal as is seen experimentally£2i21i. The source of the incoherence stems from the self-energy shown 
in Fig. (UI) which remains constant at the Fermi level at n — 0.97. This leads necessarily to a violation of Luttinger's 
theorem. In fact, the Fermi surface (defined by the maximum in the spectral function) volume at n = 0.97 is roughly 
30% larger than the Luttinger volume. In the overdoped regime, the self-energy has the characteristic dependence 
of a Fermi liquid and hence we recover Luttinger's theorem as the sharp spectral features in Fig. H17|l reveal for 
n = 0.30. Our results indicate a smooth crossover between the lightly doped regime and overdoped regimes where 
Luttinger's theorem is reinstated. 
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FIG. 18: Fermi surface in the static approximation for four fillings as indicated and three different values oiU: 1) solid line, 
U = 8t, 2) dashed line, U = lOOOt and 3) dashed-dotted line, f/ = 0, the non-interacting limit. 



However, broad spectral features are not the only contributor to the violation of Luttinger's theorem. Consider the 
static approximation in which the self energy in Eq. is explicitly set to zero. The details of this level of theory are 
derived in Appendix E. At this level of theory, the spectral function for the LHB and UHBs correspond to a series of 
(5— functions. Nonetheless, the bands generated do not describe Fermi liquid quasi-particles because each k— state still 
has spectral weight both below and above the Fermi level. Relative to the dynamical results, we find that the topology 
and volume of the Fermi surface do not change as revealed by Fig. p8|) . The solid line corresponds to ?7 = 8i, dashed 
line to J7 = 1000< and dashed-dotted to [/ = 0. Clearly shown in Fig. p8|l is the evolution from a hole to an electron- 
like Fermi surface at critical doping levels of .791 for U = 8t and .668 for U = lOOOi. The critical concentrations at 
which the curvature of the Fermi surface changes sign corresponds to Xc — 0.668 and 0.791 for U = lOOOt and U = 8t, 
respectively. As anticipated, these values are less than Xc — 2/3. While for U = 8t, Xc is remarkably close to the 
a^crit = .19 of Loram*"^, it is unclear whether the closing of the pseudogap is always accompanied with a sign change 
of the Hall coefficient. However, a sudden sign change would certainly explain the appearance of the peak in the 
density of states shown in Fig. (|10|l once the pseudogap vanishes. Three additional features are apparent. First, at 
small concentrations, regardless of U, all Fermi surfaces are electron-like and coincide with the non-interacting limit. 
An analytical proof of this result in given in Appendix E. Second, at intermediate fillings, the Fermi surface in the 
interacting system is hole-like as opposed to electron-like in the non-interacting system. Finally, the area of the FS 
for U — 8t and n — 1 — x — .97 is clearly larger than that dictated by Luttinger's theorem, 27r^(l — x) — 1.947r^. 
From the maximum in the spectral function, we find that the experimental value for the FS area in LSCO^^ for 
n — .97 is 2.067r^ which represents a non-trivial 8% deviation from the Luttinger result. Such a large deviation cannot 
be attributed to experimental uncertainty (at most 2-3%)^. Hence, while the FS computed here is clearly larger 
(2.57r^) than the experimental value, both are in qualitative agreement that Luttinger's theorem is violated in the 
underdoped regime. Of course, improved quantitative agreement with experiment can be obtained by including band 
parameters such as a next-nearest neighbour hopping interaction, t' and a doping dependent U. Hence, we see that 
even though the UHB and LHBs are sharp, the Luttinger volume is not preserved. In Appendix E, we show explicitly 
that the source of this breakdown stems from the bifurcation of the spectral weight of each k— state into a high and 
low energy part. Consequently, removing a single electron is no longer accomplished simply by removing a single 
k— state as a result of the breakdown of the band insulator sum rule. Hence, a key consequence of Mottness in 2D is a 
violation of Luttinger's theorem for n 7^ 1 as additional extensive numerical work attests^^i^^i^S42ii2S. In the heavily 
overdoped regime (see Fig. H18|) . the spectral lies predominantly in the LHB and hence one hole=one k— state and 
Luttinger's theorem is reinstated. Consequently, under hole doping, the hole and electron regimes are fundamentally 
asymmetrical as emphasized by Hirschi2i. 
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FIG. 19: Hall coefficient as a function of doping using the static approximation for the Hubbard operators with U = lOOOt. 
The inset shows that there is an antisymmetry between electron and hole doping. In the inset, the solid line corresponds to 
U = lOOOt and the dashed line to i7 = 0. Both show clearly that the deviation from 1/x is induced by the sign change rather 
than a liberation of charge. 



To compute the Hall coefficient, 



Rh ^ <7xvz/<jIx (64) 



,102 



we work within Boltzmann transport theory in whicl 



XIJZ 



^^^E--(^xV^V.(-^) (65) 



k 



(66) 



Here, 1/r is the scattering rate, 17 the volume, and / the Fermi distribution function. Our use of the Boltzmann 
equation should suffice as long as the interactions dominate which is certainly true in the case of interest, T <ti U. In 
using the Boltzmann approach, it is easiest to work in the large U limit because in this case, the spectral function 
is independent of momentum. Consequently, we used the static energy bands for U = lOOOt and compute Rh using 
Eqs. H64|l - H66|) . Fig. H19|) demonstrates that the sign of the Hall constant is consistent with the curvature of the Fermi 
surfaces shown in Fig. (|18|l . In addition, the deviation from 1/x in the region close to Xc = 1 — .668 is tied to the 
impending sign change. As the inset illustrates, Rh diverges at half-filling and changes sign for both electron and 
hole doping in contrast to weakly-interacting scenarios which can yield at most one sign change (dashed line) and no 
divergence at n = 1. Of course, the static approximation does not include the pseudogap found earlier. Note that 
regardless of which model is used for the pseudogap, the (tt, 0) regions of the Fermi surface become gapped. Unless 
the curvature of the Fermi surface is modified by the removal of the (7r,0) regions, the pseudogap cannot change the 
sign of the Hall coefficient nor eliminate the divergence at half-filling. In fact, at the doping level {x « Xc) at which 
the removal of the (tt, 0) regions is most likely to affect the curvature of the FS, the pseudogap vanishes. 



IV. FINAL REMARKS 



We have explored here a dynamical method which incorporates the local physics of a doped Mott insulator. The 
summary of our findings is catalogued in Fig. (|20|) Physics on the Mott gap scale U as well as the nearest-neighbour 
interaction scale, J, play several key roles. The Mott scale, U, sets the energy range for spectral weight transfer 
and leads to a breakdown of the band insulator sum rule. This ultimately leads to a Fermi surface volume that 
exceeds that dictated by Luttinger's theorem in the underdoped regime. Reinstatement of the Luttinger volume in 
the overdoped regime points to a fundamental asymmetry between holes and electrons in a hole-doped Mott insulator. 
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FIG. 20: Heirarchy of energy scales and the corresponding physical processes that obtain in a doped Mott insulator. AF 
represents antiferromagentic order. 



An additional role played by the Mott scale is the generation of a hierarchy of interactions of increasing range. The 
most important of these is the nearest-neighbour interaction, J « /U . Spectral weight transfer across the Mott 
gap points to an inseparability of high and low energy scales. Hence, it is unclear in what sense a true low-energy 
theory can be formulated for a doped Mott insulator. We propose that the antiferromagnet that forms in a Mott 
insulator is distinct from a spin-density wave antiferromagnet. Finally, we have found that the J is also responsible 
for the pseudogap. The pseudogap simply reflects the restricted phase space that strongly correlated excitations on 
neighbouring sites encounter. The current pattern shown in Fig. H15(l arises from such neighbouring correlations and 
could explain the origin of the direction-dependent t-reversal symmetry breaking observed in the normal states of 
the cuprates. Three-site correlations, which lead to a doping dependent spin exchange interaction, are crucial to the 
vanishing of the pseudogap at Xcrit- Because this state of affairs obtains beyond ccopt, our proposal resonantes with 
that of Loram and coUeagues^^. Finally, our work suggests that doping a Mott insulator gives rise to a heirarchy of 
energy scales all derived from the Mott gap, U . Hence, the Mott state found here has the high energy scale needed to 
explain the spectral weight transfer from 2eV to the Fermi energy when superconductivity obtains^^i^^i^. Whether 
the emergence of successively lower energy scales as a function of doping can be formulated within a renormalization 
group scheme remains an open question in strongly correlated electron physics. Nonetheless, it is along these lines 
that our current work is directed. 



V. APPENDIX A: EQUATIONS OF MOTION FOR THE TWO-LEVEL OPERATORS 

Consider the orthogonalized basis for the two-site problem, $„, Eqs. H17I18|1 . The equations of motion 

i^$„ = [$„,i7], (67) 

where H is the Hubbard Hamiltonian, can be easily obtained using the equations of motion for the single-site level 
operators, Eq. (|26|l . Assuming that <&„ is defined on x and x' , a superscript a will indicate a sum over all the 
nearest-neighbors of x with the exception of x' divided by the total number of nearest-neighbors. Similarly, a' will 



27 



include all the nearest-neighbors of x' with the exception of x. Explicitly, we have: 

- y (c!.^ + c^^'^)DBs - ficl% - c%)DBa, (68) 

^l-FBl = (2.o-.+ ^)i^i^l-^(c«-cr)Ei.+ -;=(ct<^-cf)FF^ 
+ \{c% + cl%)FFs - lidl - cI^:)FFa 

- ^{c% - cli)DBs - ^{cll + cli)DBA, (69) 

i^^FD^s = (2^0 - 3m + C/ + ^) - ^{cl% + cli)DD - ^{cl% - cli)FF- 
+ ^{c%-c^,)FFs-^{c% + c%)FFa 

- \{c% + c^,)DBs + \{c% - c^„)DBa, (70) 



'it^^^ = (;^^o-^l^ + U-^^yDl-^{c^^„-c^^i)DD + ^{c'Z+c%)FF'' 



- ^{c% + c%)FFs + ^{c%-c%)FFa 

- 1{c1,-c%)DBs + 1{c% + c%)DBa, (71) 

i^^BB = 2eoBB - 1= Y.^cf + 4" )FB5 - ^(4« - cf)FB\, (72) 
d 

i-Q^DD = 2{eo -2ij, + U)DD 

i^^FF-^ = 2{e^-,,)FF-^-l={ct-4)FB% + l={4 + 4)FBl 

- ^{cl% - c%)FD% + ^{cll + cl%)FD% (74) 

IqIFFs = 2{eo-ti)FFs 

- \ Y.^c%-c%)FBl + \Y.^c^r+c%)FB\ 

T T 
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i^FFA = 2{eo - tJ-)FFA + ^DBs 



dt 



T T 

^ ^(ct'^ + cf )FDl + ^(4" - 4">^A. (76) 



i^^DBs = {2eo -2fi + U)DBs+*-FFA 

r r 



= {2so-2^i + U)DBa 

^(^- - ^'^r)FB^s - ^ E ^(^-- + c^J^^Sa 

r T 

+ |5](4''-4")^^^5 + |E('^^'' + '^^")^^^- (78) 



VI. APPENDIX B: SELF-ENERGIES FOR THE TWO-SITE RESOLVENTS 

The resolvents associated with the two-site operators $„ can be expressed in the form 

i?$„$„(ti') = {lu — — E$„$„^(ti'))^"^ . (79) 

For the diagonal components, we use the notation = X^^ , where X is R, E, or S. We evaluate the self-energies 

S$^$,^(tj) within a one loop approximation. Using the spectral functions p+ and , Eq. of the irreducible 

propagators , we obtain 

P f 

T,FBsi^) = T / dx[2p+{x)il-f{x))RBBi(^-x)-p+{x)f{x)RFFADBsi'^ + x) 



4 „ 

+ p+ {x)f{x)RDBs {i^ + x) + p_ {x)J{x)Rdba ('^ + 

+ 3p-{x)f{x)RFFs{(^ + x)+p+{x)f{x)RFFA{'^ + x)], (80) 



5^J^S^(w) = — j dx [2p^{x){l - f{x))RBB{i^ - x) + p^{x)f{x)RFFADBs{'^ + x) 
+ P-(x)f{x)RDBs{<^ + x) + P+{x)fix)RDBAi^ + 

+ 3p+ix)f{x)RFFsi^ + x)+P-ix)f{x)RFFAi^ + x)], (81) 



^FDsi^) = ^ / [2/9+(a;)/(x)i?£)£)(tJ + x) + p+(x)(l - /(a;))i?FFAnss(^ - 2^) 

+ P-K(a;)(l - f{x))RDBs{^ - x)+ P-{x){l - f{x))RDBA{^ - x) 

+ 5p-{x){l~fix))RFFsi^-x)+p+ix)il-fix))RFFA{^-x)], (82) 



29 



P f 

SfD^(w) = I dx[2p^{x)f{x)RDD{<JJ + x)- p-{x){l- f{x))RFFADBs{'^-x) 

+ P-{x){l - f{x))RDBs{<^ - x) +P+(a;)(l - f{x))RDBAi!^ - x) 

+ Zp+{x){l-f{x))RFFs{^-x)+p-{x){l-f{x))RFFA^-x)], (83) 
^bb{co) = P j dx [p+{x)f{x)RFBsi^ + x)+P-{x)f{x)RFBA{'-J + x)], (84) 



1:dd{oj)^P J dx [p+{x){l - f{x))RFDsi(^ - x) + p-{x){l - f{x))RFDA'-^ - x)], (85) 

P r 

T.FFsi'^) = y / dx[p-{x)f{x)RFDs{'^ + x)+p+{x)f{x)RFDA{^ + x) 
+ p_{x){l-f{x))RFBs{'^-x) 

+ p+{x){l- f{x))RFBA{^ + x% (86) 

^FF' (W) = ^FFs {^) = ^DBa (w) (87) 

P r 

^FFa{^) = / dx[p+{x)f{x)RFDs{^ + x)+P-{x)f{x)RFDA{'^ + x) 
+ P+{x){'i-- f{x))RFBs{(^-x) 

+ P-{x){l-f{x))RFBA{'^ + x)], (88) 



p 



T^FFADBsi^) = 2" y dx [p+{x)f(x)RFDs{^ + x) - P-{x)f{x)RFDA{'^ + x) 

- P+{x){l- !{x))Rfbs{^-x) 

+ p_{x){l- f{x))RFBA{^ + x)\, 



^FFaDBs{^) = '^DBsFFa{^), (90) 

S,3SsM = SfF^M- (91) 

VII. APPENDIX C: SELF-ENERGIES CORRECTIONS DUE TO SPIN FLUCTUATION 

We include the effects of spin fluctuation as higher order corrections to the self-energies of the resolvents. The 
spectral functions associated with local singlet and triplet states are sharply peaked at well defined energies separated 
by the effective antiferromagnetic coupling constant 



/oo 
-00 



o-FFs - ^FFa) (92) 
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The effects of spin fluctuations with the environment can be approximately described by the eff'ective antiferromagnetic 
interaction 

5H^ff = ^Jnn«. (93) 

where n = cVc is the spin density, cr, are the Pauli matrices and = (c|,c|). Although a proper account of the 
singlet-triplet mixing can be given by considering vertex corrections, we can use a simpler two-step approach: 1) 
Working in the basis formed by the eigenstates of the two-site problem (as we, in fact, do) we take care of the spin 
fluctuations for the cluster without a bath. As we introduce the bath, an additional singlet-triplet mixing occurs 
which is not captured by NCA (as shown by the very sharp features in the FFA and the FFS resolvents). We try 
to approximate this mixing by an effective spin-spin interaction. As a result of this mixing, the FFA and FFS states 
are broadened. 2) Due to the the self-consistency of the approach, this spin-spin interaction is also present between 
the cluster and the bath. The corrections that we introduce are the effect of this additional effective interaction with 
the bath on the self-energies of the resolvents. Consequently, some of the equations of motion for the two-site level 
operators (see Appendix A), will be modified. Explicitly we have: 

i^FB% = ^Knf +nf')FB§ + a(nf-nf')FB5 

+ {n% + n%)FBs^ + in%-n%)FB^''] + ---, (94) 

i§-^FB^A = ^[a{nl-nf)FB^, + a{nl+nf)FBl 

+ {n%-nOFBs^ + {n%+n^_:,)FB^'^] + ..., (95) 

'F/^^' = ^[a(nf +nf)FZ?g + a(nf-nf)FZ?^ 

+ {n%+nOFDs^ + {n%-nOFD^-] + --- (96) 

'§iFD^A = ^[a{nl-nf)FD^s + a{n^+nf)FD'^^ 

+ {n%-n'l,)FDs" + {n%+n%)FD^''] + ..., (97) 

i^FF- = J[a(nf + nl')FF^ + ^{n% + n^FFs 

- ^{n%-nOFFA] + ---, (98) 

i^^FFs = J[a{n1-4)FFA + ^{n% + n'^,)FF'^ 

+ -L(n«,+n^;)Fi^-'^] + ..., (99) 

i^FFa = J[ainl-nf)FFs-^{n%-n^_:,)FF^ 

+ l={n%-n''_:,)FF-'^]+--., (100) 

where we used the notation =ni+ ian2- These terms will determine the self-energy corrections 

SEpBsi^) = J dx ^J'^[f+{x){l ~ f{x))RFBs{^ - 3:) 

+ ^.{x){l-f{x))RFBA^-^)], (101) 
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SY.pbA'^) = j dx^J^[ip_{x){l- f{x))RFBs{^-x) 

+ ^+{x){l- f{x))RFBA^-x)], (102) 



SY.FDsi'^) = j dx'^j'^[(p+{x){l - f{x))RFDs{'^ - 

+ ip.{x){l- f{x))RFDA^-x)\, (103) 



SIIfdA'^) = j dx^J^[ip_{x){l- f{x))RFDs{^-x) 

+ ^+{x){l- f{x))RFDA^-x)], (104) 

SllFFsi'^) = j dx 2J^[2ip+{x){l - fix))RFFA^ - x) 

+ <p_{x){l-fix))RFFA^-x)], (105) 



ST.ffA'^) ^ J '^^ 6JV-(a;)(l - I{x))Rffs{^ - x), (106) 
where ip±{y>) = fiui) ± ip'{to) and 

= F.T.{{nl,nl}) = ^F.T.{{n%,n%}) 

^'{oj) = F.T.{{nt.n'i}) = \F.T.{{n^,,n%}) (107) 

Note because <p(a;) and </''(w) are expressed in terms of functions defined on a and a' which denote the nearest 
neighbours of the two sites in the cluster, we have avoided double counting the spin fluctuations on the cluster. The 
functions ip{Lo) and i^'i^o) arc computed within a non-crossing approximation, taking into account the only dominant 
contributions (containing products of on-site and nearest-neighbor operators). Introducing the coordination factors 
A = {2d-l)/{{2df) and A' = (2d - 2)/((2d)2), we obtain 

ip{u)) W ^ J ^^{\ t^^^® {x)^FBs {U}-X)+ GFBs {x)(JFBs (<^ " x)] 
+ \ V'FBs {x)aFBA {l^- X)+ OFBs {x)(TfBa (^ - x)] 

+ ^ Wfba {x)^fbs (w - a;) + aFBA {x)(Tfbs - x)] 
+ \ Wfba {x)^fba (w - a;) + gfba {x)crFBA - x)] 

+ ^ [f^FDs {x)aFDs {ui - X)+ ffFBs {x)(TfBs ('^ " x)\ 

+ ^ [(^FDs {x)^FDa {cO- x)+ a-FBs {x)crFBA " x)] 

+ ^ WfDa {x)&FDs (w - x) + ^fBa {x)<TFBs " x)] 

+ ^ WfDa {x)^FDa {^^ - x)+ aFBA {x)<JfBa (^ - x)] 
+ 2[aFFs{x)aFFs{'^ - x) + aFFs{x)(JFFs{'^ " a;)] 

+ 0-FFs{x)aFFA{^ - x) + aFFsix)<yFFAi^ - x) 

+ C^FFAix)aFFsi'^ - x) + aFFA{x)(^FFs{'^ - x)^ , (108) 
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A' 
Z 




apBs {iJJ - x) + (TFBs (x)' 


<7FB, 


1 
2 


V'FBs {x)aFBA ('^ - 


x) 


- d-FBsix)aFBA - 


x)] 


1 
2 


[o-fba {x)^fbs - 


x) 


- O-FBa {x)aFBs - 


x)] 


1 

+ 2 


[<^FBa {x)aFBA - 


x) 


+ ^FBa{x)(TfBa{'^ - 


-x)] 


1 

+ 2 


[(JFDs {x)dFDs - 


x) 


+ &FBs{x)(^FBs - 


x)] 


1 

^ 2 


[crFDs{x)^FDA - 


x) 


- 0-FBsix)crFBA - 


x)] 


1 
2 


WfDa {x)&FDs ('^ - 


x) 


- o-fba {x)aFBs - 


x)] 


1 

+ 2 


WfDa {x)^FDa - 


x) 


+ ^FBa {x)(TfBa{'^ - 


-x)] 



+ 2[(JfFs{x)(^FFs{^ — x) + aFFsix)(JFFs{^ — x)] 

- aFFs{x)^FFA{^ " x) ~ aFFs{x)cTFFA{^ ~ x) 

- aFFA{x)^FFsi^ ~ x) - aFFA{x)<yFFsi^ - x)^ ■ (109) 

Note that in these equations, it is the Fermi function which appears not the Bose distribution function. In general, the 
expression for the self energy depends on which form one uses for the spectral function. We used F.T.l{n3{t),n3{t')}{= 
ip{uj). Usually the pectral function is defined by F.T.{[n3{t),n3{t')]) = p{uj). They are related as (p{uj) = [{e^'^ + 
\)l{e^'^ - The loop contains 



which, in terms of ip^ becomes 



g/3(L^+C^') 

du}' ip{J ^ u:)R{;^J) ^ -, (111) 



thereby justifying the Fermi function in the self-energy corrections. 

VIII. APPENDIX D:LEVEL OPERATOR REPRESENTATION OF 5 J AND EXPLICIT EXPRESSIONS 

FOR Dmo AND Dmi 

With the assumption used to derive Eq. p2fl . we find that the the dynamical correction operator 5J , Eq. JTJ), can 
be expressed in terms of two-site level operators, 

(112) 





5Ja 


= ~i 










With the notation 














ai 




















n 

6. = -H 








we obtain for the coefficients 














flSB FB%, — 


-1. 


-6i), 


'^FD^' DD ~ 






-&2): 


aBB FB". — 

A 


-1. 


+ foi), 


'^FD^" DD ~ 




+ a2 


-62). 


^FB-" FFs ^ 


1, 

2(«i + 


bi), 


flFFs FD^ = 






-^2), 


^FB^" FFs ~ 


-1. 


-bi), 


OFFs fd^ — 




- 02 


-?>2) 
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O-FB^ 


FF" 




Clpp-a FDg" ~ 




+ 02 


- &2), 




FF" 




O-FF-" FD^" = 


5" 


- 02 


-b2), 




FFa 


= cr(l + fli - bi), 


a_FFA -FDg = 


id. 


02 + 


b2), 


"■FB^" 


FFa 


= a{l - ai - 61), 


O-FFa FD'^ = 


^<'- 


02 + 


b2), 


aoBs 


FD% 




^FBg" DBs ^ 




- 02 4- 


-b2), 


aOBs 


FDl 




'^FB^" DBs ^ 




- 02 


+ b2), 


O-DBa 


FDl 




^FB^' DBa ^ 




02 - 


62), 


O-DBa 


FDl 


= ^(ai-&i), 


^FB^' DBa ^ 


— cr , 


- 02 


-62). 
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Using these coefficients, we can write explicitly the dynamical corrections DrriQ and Dmi. From Eq. (|12|) . we have 
Dmo(w) = % I dxdx' ] , , (115) 

^ FBs\-'^FBsi^)^BBix') + d-FBs{x)'^BB{x')] 

+ 0?BB FBa\-'^FBa{x)^Bb{x') + aFBA{x)<yBB{x')\ 

+ o^FDs ddWdd{x)&fDs{x') +^dd{x)(tfDs{x')] 

2 



+ Oj?£,^ dd\'^dd{x)cffda{x)'^odd{x)ufda{x)\ 

+ ("fBs -FFs+'^FSs FF'')[^-Fi='s (2^)^-FSs (2;') + ^F_Fs(2^)^-FSs (2;')] 

+ WfBa FFs +0-FBA FF'')['^FFs{x)&FBa{x') + &FFs{x)(TFBa{x')] 

+ (offs fds + Off- fds ) I'^^'C's (2;)o-ffs {x) + ^fds (a;)crFFs (x')] 

+ (ofFs FDa + «FF- FDa ) I'^^'Oa (2;)o-FFs (^^O + ^FDa {x)(TfFs {x')] 

+ a],Bs FFa \-^FFa (xpFBs (x') + ^FFa {x)(JfBs {x')] 

+ O^FBa FFa \^FFa {x)ofBa {x') + ^FFa {x)(JFBa {x')] 

+ '^FFa FDsl'^FDsix)o-FFAix') + 0-FDsix)o-FFAix')] 

+ a%p^ WfDa {x)^fFa {x') + ^fDa {x)(Tfpa {x')\ 



+ 2apBs DBs'^FBs FFaWdBs FFa{x)&FBs{x') +<^DBs FFa{x)(JfBs{x')] 

+ '^O-FBa DBsO-FBa FFaI'^DBs FFa{x)^FBa{x') + ^dBs FFa (2;)o'FBa (2^')] 

+ 2aDBs FDs'^FFa FDs['^FDs{x)aDBs FFa{x') + &FDs{x)(^DBs FFa{x')] 

+ 2aDBs FDaO'FFa FDa[^FDa{x)&DBs FFa{x') + &FDa{x)(^DBs FFa{x')\ 

+ "-FBs DBs i^DBs {x)&fBs [x') + &DBs {x)<7pBs {x')] 

+ OfBa DBs ['^DBs {x)^FBa {x) + &DBs ix)crFBA (x')] 

+ O-DBs FDs [^FDs {x)^DBs {x') + &FDs {x)(JdBs {x')] 

+ O-BBs FDa I'^FDa {x)^DBs {x') + ^FDa {x)(TDBs {x')] 

+ O^FBs DBa \-^DBa {x)^FBs {x') + ^DBa {x)(JpBs {x')] 

+ a],BA DBa WdBa ix)^FBA (x') + ^DBa (x^FBa (x')] 

+ O-DBa FDs i'^FDs {x)&DBa {x') + apDs {x)<7dBa {x')] 

+ OdBa FDa ['^FDa {x)^DBa {x') + ^FDa {x)<JdBa {x')] ] ■ 
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Similarly, for Dmi we obtain 

P r 1 

Dmi(uj) = — / dxdx' r (116) 

Z J uj-x + x' + tS ^ ^ 

X {flSS FBsl-'^FBs{x)aBB{x') + aFBs{x)<JBB{x')] 

- O.BB FBAl-'^PBA{x)d-BB{x') + d-FBAi^^BBix')] 
+ O-FDs DD['^DD{x)dFDs{x) + ^Dd{x)(TFDs{x)] 

- '^FDa DD['^DD{x)aFDA{x') + &DD{x)aFDA{x')] 

- (o-FBs FFs + O'FBs FF" ) [^FFs {x)^FBs (x') + ^FFs {x)<JfBs i^')] 
+ (apBA FFs + ^-FBa FF'' ) [f^FFs {x)^FBa i^') + ^FFs {x)<TfBa i^')] 

- {a^FFs FDs +(^'fF-' FDs)\-^FDs{x)&FFs{x') + aFDs{x)(TFFs{x')\ 
+ {apFs FDa +"'FF' FDA)\-'^FDA{x)aFFs{x') + &FDa{x)<^FFs{x')\ 

+ apg^ FFA\-^FFA{x)aFBs{x') + &FFa{x)(^FBs{x')] 

- O.FBA FFa i'^FFA {x)^FBa i^') + ^FFa {x)(JfBa {x')] 
+ O^FFa FDs i^FDs {x)&FFa {x') + ^FDs {x)<yFFA {x')] 

- O^FFa FDa I'^FDa {x^fFa (x') + &FDa {x)(JFFa {x')] 

+ 2aFBs DBs'^FBs FFaWdBs FFa{x)&FBs{x') + auBs FFa{x)(JfBs{x')\ 

- '^O.fBa DBsO-FBa FFa[<^DBs F Fa{x)^ F Ba{x') +^DBs FFa (2;)o'FBa (^^O] 
+ '^a.DBs FDs^'FFa FDs[^FDsix)o'DBs F Fa{x') + ^ F Ds{x)cF DBs FFa{x')] 

- '^aoBs FDaO'FFa FDa[^FDa{x)^DBs FFa{x') + ^FDa{x)(JdBs FFa{x')] 
+ (^FBs DBs i^DBs {x)&FBs {x') + &DBs {x)<7FBs (x)] 

- jjg^ [aoBs {x)^fBa (x) + ^dBs {x)<JfBa (x')] 

+ O.DBS FDs [^FDs {x)^DBs {x') + &FDs {x)(JdBs {x')] 

- OdiBs FDa [^FDa {x)^DBs {x') + ^FDa {x)(TdBs {x')] 

- O^FBs DBa [<^DBa {x)&FBs {x') + ^DBa {x)<7FBs (x')] 
+ a,FBA DBA WdBa {x)aFBA {x') + aoBA {x)(TfBa {x')] 

- ai)i3,4 FDs V'FDs {x)^DBa (x') + ^FDs {x)(JDBa (x')] 

+ -aoBA FDa i'^FDA {x)aDBA (x') + ^FDa {x)(yDBA (x')] } ■ 

IX. APPENDIX E: STATIC APPROXIMATION 

We show explicitly in this appendix that the Hubbard operator technique in the static approximation^'^'^ correctly 
recovers the non-interacting or Fermi liquid limit when U ~ regardless of the filling and as n — > for any U. 
Consequently, the violation of Luttinger's theorem found here is not an artificat of the method. Consider the Hubbard 
operator basis, ipi = ^ia- and -02 = Via- and the associated Green functions Sap = {{'ipa]''Pi3))- Within the static 
approximation^"'^, the expression for the retarded Green function in Fourier space becomes 

2 cr^^^ (k) 
'^-^ uj — e,(k) + 10 

The dispersion relations for the two bands are given by ei,2(k) = i?(k) ± Q(k) with 

i?(k) = lu-ti--^[m(k)+8ta{k)hh], (118) 
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■h)] 



where Ji = 1 — n/2 and I2 = n/2 and Q!(k) = ^(cos(fca;) + cos(A:y)) and we used the notation TO(k) = 4t[e + Q!(k)(p- 

and5(k) = -f/ + (l-n)/(/i/2m(k)) ^he^e e = {^Uf,) - {r^M and p = (ni.<J + (5,5|")- (6^61") with^i 
bi = Ci-[Cii and a indicates a sum over nearest neighbours of site i. The explicit expressions for the spectral functions 

^ali^) are given by 



(k) = ^ 



1 + 



5(k) 



(k) 



(k) 



m(k) 
2Q(k)' 

2 



1 



2Q(k) 



5(k) 



5(k) 



2Q(k), 



J2) 



(k) 



2Q(k) 



m(k) 

2g(k) 

h 

1 + 



(120) 



5(k) 



2<3(k), 



Note the fact that the spectral functions are k-dependent and they also depend on the doping level, through Ii and 
I2, and temperature, due to the self-consistent parameters e and p. At half- filling, the spectral function for the lower 
Hubbard band determines the amplitudes in |MI) through = {uk + Vk)^/2. In the strong-coupling limit, U t, 



the only dependence that remains is on filling. In this limit a^^^ 
All of our calculations of the Hall coefficient were performed in this limit. 

Consider now the two weak-coupling limits in which Fermi liquid theory should hold: Case a) n ^ and Case b) 
f/ — > 0. In the first case, g — > —U + 2m/n, Q — > 131/2, the correlations in p become independent and hence p oc v?, 
e cc n implying that m(k) cx n. Consequently, the dispersion for the lower Hubbard band reduces exactly to that of 
the non-interacting limit, ei(k) = — — 4to(k) = eo(k). Moreover, all the spectral weight resides in this band because 
as n — > 0, g/2Q = 1 and m(k) — > 0, implying that a^^^ = 1 and a^^^ = 0. That the static approximation reduces to the 
correct non-interacting limit is not unexpected as Figures (1) and (2) illustrate that the Fermi surface is independent 
of J7 as n ^ 0. In the U ^ limit, (/(k) (1 — n)m(k)//i/2 and as a consequence, Q = |m(k)|/2/i/2. As a result, the 
band dispersion relations are ei,2 = eo(k) — (m(k) =f |m(k)[)/2Ji J2 with spectral weights cr^^'^^ = l/2±m(k)/(2|m(k)|) 
which are either unity or zero. Consequently, although two bands still exist, only the free particle dispersion carries 
unit spectral weight because the |TO(k)| terms enter with opposite signs. Hence, the static approximation correctly 
reproduces the non-interacting limit when U 0. As a consequence, the violation of Luttinger's theorem seen here 
is not an artifact of the approximation scheme but stems fundamentally from the splitting of the spectral weight over 
two bands although no symmetries are broken, the hallmark of Mottness. 
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